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ABSTRACT 

We present the derivation of the free-free emission on the Galactic plane between 
£ = 20° and 44° and |5| < 4°, using Radio Recombination Line (RRL) data from 
the Hi Parkcs All-Sky Survey (HIPASS). Following an upgrade on the RRL data 
reduction technique, which improves significantly the quality of th e final RRL spectra , 
we have extended the analysis to three times the area covered in lAlves et al.l ()201Cll ). 
The final RRL map has an angular resolution of 14.8 arcmin and a velocity resolution 
of 20 km s"^ 

A map of the electron temperature (Tc) of the ionised gas is derived for the area 
under study using the line and continuum data from the present survey. The mean Tc 
on the Galactic plane is 6000 K. The first direct measure of the free- free emission is 
obtained based on the derived Tg map. Subtraction of this thermal component from 
the total continuum leaves the first direct measure of the synchrotron emission at 
1.4 GHz. A narrow component of width 2° is identified in the latitude distribution of 
the synchrotron emission. We present a list of Hii regions and SNRs extracted from 
the present free-free and synchrotron maps, where we confirm the synchrotron nature 
of three objects: G41. 12-0.21, G41. 15-^0.39 and G35. 59-0.44. We also identify a bright 
(42 Jy) new double radio galaxy, J1841-0152, previously unrecognised owing to the 
high optical extinction in the region. 

The latitude distribution for the RRL-derived free-free emission shows that the 
WMAP Maximum Entropy Method (MEM) is too high by ~ 50 per cent, in agreement 
with other recent results. The extension of this study to the inner Galaxy region 
£ = —50° to 50° will allow a better overall comparison of the RRL result with WMAP. 

Key words: radiation mechanisms: general - methods: data analysis - dust, extinc- 
tion - Hii regions - ISM: lines and bands - Galaxy structure - radio lines: ISM 



1 INTRODUCTION 

Radio Recombination Lines (RRLs) are a powerful tool 
for the diagnostics of the interstellar medium, tracing the 
ionised component, its electron temperature, velocity and 
radial distributions. They can be used to determine the 
free-free Emission Measure on the Galactic plane where 
heavy obscuration makes optical measurements impossible. 

* E-mail:marta. alves@ias.u-psud.fr 



In lAlves et al] l|2010l ). hereafter referred to as Paper l, we 
have presented the unique method of determining the free- 
free emission on the Galactic plane using RRLs. In the 
present work we include a determination of the electron tem- 
perature of the ionised gas, following the separation of the 
continuum free-free and synchrotron emission. 

Once the free-free emission is determined, the syn- 
chrotron component is available by comparing with the total 
continuum emission at the same frequency, from the same 
survey. This approach is more rigorous than the normal 



© RAS 



2 M.I.R. Alves et al. 



method of using the difFerent spectral indices of the two 
components to make the separation. It allows a first clear 
image of the Galactic latitude and longitude distributions of 
the two components. 

Paper I showed the feasibility of the RRL approach. 
The aim of the present paper is to extend the work to a 
larger longitude range in order to sample a greater volume of 
the Galaxy. The original paper covered the Sagittarius arm 
in the first quadrant; here we include more of the Sagit- 
tarius arm and a Section of the Scutum arm. Separating 
the emission from the two spiral arms allows the investi- 
gation of the z-thickness of the ionised layer. Early deter- 
minations of the latitude distribution of the diffuse RRL 
emission on the G alactic plaiic gave a thi n layer of scale 
height 40 - 100 pc iGordon fc Catol (|l97l ): iHart fc Pedl"a3 
l|l976t ): [Meige3 (f IQTsI ). These values are in agreement with 
the z-distribution of the Hll regions within the solar cir- 
cle, with a F WHM (Full Width at Ha lf Maximum) of 
90-120 pc l|Paladini. Davies fc DeZott i 2004) and some- 
what greater than that of t he OB stars with a FWHM of 
32 pc ({Bronfman et al.ll2000l ). 

The analysis of Paper I was based on the deep ZOA 
Galactic plane survey which made a first attempt at re- 
moving the spectral ripples in RRL work. We improve the 
signal-to-noise of the spectra and the spatial ripple by in- 
cluding the HIPASS data which was acquired by scanning 
in declination, as opposed to Galactic longitude scanning in 
the ZOA survey. We describe the update of the RRL reduc- 
tion pipeline used in Paper I, which includes a new bandpass 
algorithm exclusively designed for the analysis of RRL spec- 
tral line data in HIPASS. It extends the analysis from the 
previous 8°(^) x 8° (6) region centred at {i,b) = (40°, 0°) to 
a 24° (^) X 8°(fe) area centred at {l,b) = (32°, 0°). 

The paper is organised as follows: Section [5] gives the 
HIPASS/ZOA survey parameters. In Section [3] we describe 
the RRL data reduction and the properties of the final data 
cube. The derivation of the free-free map is presented in 
Section |4l along with a catalogue of Hll regions extracted 
from the same map. In Section [S] we present the synchrotron 
map derived from this survey with its catalogue of SNRs. 
Section [6] compares the RRL estimates of the free-free with 
WMAP data. Finally, Section [7] gives the velocity and radial 
distribution of the ionised gas in the Galaxy. 



2 THE RRL SURVEY 

The sensitive survey described in this work is a 
product of the Hi Par kes All-Sky Survey (HIPASS, 
IStavelev-Smith et al.l 1996f) and associated Zone of Avoid- 
ance (ZOA) survey l|Stavelev-Smith et al.lll998l ). The sur- 
veys searched for Hl-emitting galaxies in the redshift veloc- 
ities of —1280 and 12700 km s~^ using the Parkes 21cm 
multibeam receiver. Within the 64 MHz bandwidth there 
are three RRLs: H168a, H167a and H166a at 1374.601, 
1399.368 and 1424.734 MHz, respectively. 

The HIPASS survey covers the sky south of S = +25? 5. 
The scans are taken at constant right ascension, in declina- 
tion strips of 8° separated by 7arcmin. The ZOA deep sur- 
vey is scanned at constant latitude and each 8° longitude 
scan is separated by 1.4arcmin. It covers the Galactic plane 
accessible from Parkes, i = 196° - 0° - 52° and |6| < 5°. 



Fig. \T\ shows the coverage of the 8° {£) x 8° (6) region cen- 
tred at {£,b) — (40°, 0°) analysed in Paper I along with the 
24° {£) X 8° (6) region analysed in the pr esent work, overlaid 
on the radio map at 2.7 GHz from the iReich et all l{199(i ) 
survey. 

The multibeam receiver has 13 beams set on a hexago- 
nal grid with a mean observing beamwidth of 14.4 arcmin. 
The footprint of the receiver on the sky is 1?7, thus each 
scan maps an area of 8° x 1?7. The multibeam correlator 
cycle is 5s and the scan rate is l°/min. The HIPASS/ZOA 
frequency coverage is from 1362.5 to 1426.5 MHz centred 
on 1394.5 MHz and is divided into 1024 channels spaced by 
62.5 kHz, or 13.2 km s~^ for the Hi line. Therefore an 8° scan 
is about 100 integrations. The total integration time is 450 s 
per beam for the HIPASS survey and 2100 s per beam for the 
five times deeper ZOA survey. This results in a typical rms 
(root mean square) noise of 13mJy/beam and 6mJy/beam 
for the HIPASS and ZOA surveys, respectively. 



3 THE RRL DATA 
3.1 Data Reduction 

The RRL data reduction pipeline has greatly improved since 
the first analysis described in Paper I. The previously pre- 
ferred al gorithm for bandpa ss correction in the livedat43 
package I Barnes et al.ll200ll ) was minmed5, developed for 
the treatment o f extended emission from High Velocity 
Clouds (HVCs) jPutman et all 120021 '). The refinement of 
MINMED5, the minmedn algorithm, was tested in the new 
analysis. This method is the running median in a box of 
size n integrations computed along the length of the scan, 
for which the bandpass correction is the minimum of those 
medians, minmedn was found to produce spectral negatives 
and loss of emission on the Galactic plane, similarly to min- 
med5. The spectral negatives were found to arise from large- 
scale distortions of the bandpass response by very strong 
continuum sources, for which the system temperature Tsys 
increases by nearly a factor of two. For these reasons, a new 
algorithm was developed for the exclusive treatment of ex- 
tended emission from RRLs in livedata, called tsysmin„. 

tsysmin„ locates the n integrations in a scan for which 
Tsys is minimum and, for each spectral channel, takes the 
median value over these integrations as the bandpass cor- 
rection, n = 10 provides a sufficient number of integrations 
to compute a reliable median without the averaging length 
being too extended in regions where a large part of the scan 
may be occupied by extended emission. Fig.[2]shows a chan- 
nel map of the ZOA-032 datacube at V = —13.4 km s~^, 
reduced with minmedio and tsyaminu), and H167q spec- 
tra towards G35.6-f0.0. The tsysmin estimator is used be- 
cause it is much less susceptible to baseline negatives, ren- 
ders cleaner spectra with higher peak values and recovers 
more extended emission on the Galactic plane. 

After bandpass correction and calibration the spectra 
are smoothed using the Tukey filter to suppress the Gibbs 
ringing caused by the strong Galactic Hi line, which afi'ects 
the H166a RRL. The three RRLs are extracted from each 
spectrum, shifted in velocity and stacked using a weighted 



^ http:/ /www. atnf.csiro.au/computing/software/Iivedata. html 
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Galactic Longitude 



Figure 1. The coverage of the 8°{£) X 8°{b) region analys ed in Paper I (gree n) and the larger 24°(£) X 8°{b) area under study in the 
present work (red). The map is from the 2.7 GHz survey bv lReich et al] l ll99Cl) . 
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Figure 2. (Top) Maps of the ZOA-032 datacube at V = —13.4 km s^^ reduced writh minmedio (left) and tsysminio (right). (Bottom) 
The corresponding H167a spectra at {£,b) = (35?6,0?0). 



mean, where the weights are proportional to the square root 
of the bandpass response measured at each line's rest fre- 
quency. The weights are 0.35, 0.39 and 0.26 for the H168a, 
H167q and H166a, respectively. 

The 3D £ — 6— velocity cube for each 8° x 8° zone is 
generated from the weighted median of the stacked spectra 
within a cut-off radius of 6 arcmin of each pixel, using the 
package GRIDZILLA. The weight for each input spectrum is 



based on its angular distance from the pixel and its system 
temperature. The resulting data cube has 4 arcmin^ pixels 
and a LSR velocity range of ±335 km s~^. 

Further improvement of the spectra was achieved by 
combining HIPASS data with ZOA data. The HIPASS dec- 
lination scans cross the plane at an angle of ~ 30° in the lon- 
gitude range under study here and therefore yield a real min- 
imum for the bandpass estimation, as opposed to the longi- 
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Figure 3. ZOA (left) and ZOA+HIPASS spectra after DC level correction (right) at two positions; (top) G33. 1-0.1 and (bottom) 
G28. 6-2.1. The improvement after combining both datasets is visible even for spectra away from the Galactic plane. 



tude ZOA scans. We first calculate the DC level between the 
two surveys by subtracting the ZOA cube from the HIPASS 
cube, after going through the same reduction pipeline. The 
DC offset added back to the ZOA spectra is computed as 
the median of the difference cube at each Galactic latitude 
and spectral plane. After correction, both datasets are av- 
eraged giving 5 times more weight to the ZOA survey since 
it is 5 times deeper than the HIPASS survey. Finally, spec- 
tra in the combined cube are further Tukey smoothed and 
baseline removed by fitting a polynomial. Fig. shows two 
examples of spectra improved with the DC level correction 
and combination of ZOA and HIPASS data. The negative 
at V ~ 40 km s~^ in the ZOA spectrum of G33.1-0.1 is 
removed after combined with HIPASS. The improvement of 
the baselines is clear for spectra on and away from the plane, 
as Fig. 13] (d) shows. 



3.2 Final Data Cube 

The three 8° x 8° x 51 channel datacubes are combined into 
a final 24° x 8° x 51 channel datacube that covers the £-range 
20° to 44°, l&l < 4° and Flsr = ±335 km s'^ The final spec- 
tral resolution after Tukey-I- Tukey smoothing is degraded 
from 16 km s~^ to 20 km s~^. The typical noise level in each 
channel is 3.5mJy/beam or 3.0 mK. The final spatial reso- 
lution was measured on the total integrated RRL map using 
two compact Hll regions, W49 and W40. The average beam 
size is 14.8 arcmin with an uncertainty of 0.6 arcmin, which 



is consistent with the intrinsic resolution of 14.4 arcmin. It 
is also consistent with simulations, which inject false Gaus- 
sian sources into to the data prior to imaging, where, with 
weighted median gridding, the broadening of the intrinsic 
14.4 arcmin beam is negligible. The brightness temperature 
conversion of 0.84 K/Jy is obtained from the Rayleigh- Jeans 
relation for a Gaussian FWHM (Full Width at Half Max- 
imum) of 14.4 arcmin at the frequency of the H167a line. 
This was checked using well-known bright sources and has 
an uncertainty of 0.03 K/Jy. The conv ersion factor is consis- 
tent with p revious measur e ment s by IStavelev-Smith et al.l 
(2003) and iKalberla etall l|2010l) that indicate a value of 
0.80 K/Jy. 

The total RRL emission map is obtained by integrat- 
ing the spectra in the appropriate velocity range, to include 
the different contributions from the Local {V — Q km s~^), 
Sagittarius {V = 40 km s"^) and Scutum (V = 80 km s"^) 
arms. The integral varies between —20 and 150 km s~^ in 
order to account for any velocity spread in within the 
arms. The rms of the 24° x 8° RRL integrated map is 
0.05 K km s~^, and the maximum and minimum are —0.2 
and 26.7 K km s~^, respectively. 

The data are calibrated in the main beam scale, appro- 
priate for point sources, rat her than on the full-b eam which 
includes the near sidelobes ([Reich fc Re ich 1988). However 
due to the imaging process the flux for point sources is low- 
ered in the final cubes. This depends on the cut-off radius 
used in GRIDZILLA since the final value for a given pixel is 
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Figure 4. Spectrum for the Hll region W47 {£, b) = (37? 8, -0?3) 
allowing tiie H, He and C RRLs. Tlie vertical dotted lines indicate 
the rest velocity of each line. The C RRL is amplified by maser 
emission by the foreground gas and shifted to lower velocities. 



the median of all the data that are within a certain distance 
from it. Simulations show that for a cut-off radius of 6 arcmin 
the peak fluxes of point sources are underestimated by 20 
per cent. Further investigation using the datacubes reduced 
with smaller cut-off radii, for which this effect is increas- 
ingly small, show that it only affects the two most compact 
Hll regions W49 and W40 which are about 5 arcmin in size. 
The 6 arcmin cut-off radius is used because it provides an 
adequate flux scale with an uncertainty of approximately 10 
per cent for extended sources, in a low noise cube. 

The total power map from the same survey is available 
for the region under study (Calabretta et al., in prep.) and 
is shown in Fig. [T] It is reduced using ZOA and HIPASS 
data similarly to the RRL analysis, in order to recover the 
zero levels on the Galactic plane and gridded with the same 
parameters. 



3.3 Helium and Carbon RRLs 

Helium and Carbon RRLs are observed in the 
HIPASS/ZOA spectra. The expected velocity separation of 
122.1 km s~^ and 149.5 km s^^ between the hydrogen RRL 
and the helium and carbon lines, respectively, allow their 
identification. Due to the spectral resolution of 20 km s~^, 
these two RRLs are sometimes blended. When the carbon 
RRL is amplified by maser emission by foreground gas, it 
shifts towards lower velocities and the two can be separated. 
An example is shown in Fig. 3] for the Hll region W47. 

Further investigation will provide information on the 
He/C RRLs from the diffuse ionised gas in the inner Galaxy. 
This will be presented in a future work. 



4 CONVERSION TO A FREE-FREE 
TEMPLATE 

The RRL emission from a diffuse ionised hydrogen gas can 
be expressed in terms of the integral over spectral line tem- 
perature as 



where f is the frequency (kHz), Te is the electron temper- 
ature (K) and the emission measure EM is in cm~^pc. The 
corresponding continuum emission brightness temperature 
is 

Tb = 8.235 X 10~^a{n)T-°-^^iy^^^{l + 0.08)EM (2) 

where a(Tc) is a slowly var ying function of temperatur e and 
u is the frequency in GHz (jMezger fc Henderson|[l967h . The 
(1-1-0.08) term represents the additional contribution to Tb 
from helium. These equations lead to an expression for the 
ratio of the line integral to the contirmum 



6.985x10^ 



-T~ 



Tb a(Te) n(He)/(l n(H)) 

where V is in km s"^ jRohlfs fc Wilsonll2004h 



4.1 The Te of the thermal emission from 

44° 



't'GHz(3) 



20° to 



In order to convert the RRL integral into a free-free bright- 
ness temperature, a value for the electron temperature of the 
ionised gas is needed. The variation of To in discrete Hll re- 
gion s with longitude, or Galactic radius Rp, is well-known 
(Sh aver et al.lll983l : [Paladini. Davies fc DeZotti| [2004''l. It is 
caused by the greater cooling of the Hll regions closer to the 
Galacti c centre due to a higher me t al con tent in the inner 
Galaxy. iPaladini. Davies fc DeZottil l|2004l ') uses all the pub- 
lished data for 404 Hll regions with reliable Galactocentric 
distances to derive 



Tc[K] = (4166 ± 124) + (314 ± 20)i?G[kpc]. 



(4) 



It follows that the mean Tc at the Galactic radius of the Lo- 
cal arm, 8 kpc, is 6680 K and decreases to 5740 K at 5.0 kpc, 
the radius of the Scutum arm. 

The line-to continuum ratio given by equation ((3)l is one 
of the most reliable methods to determine the electron tem- 
perature of an Hll region and this survey has the advantage 
of providing both the line and continuum data at the same 
resolution, after going through the same reduction pipeline. 
However, Tb in equation (|3]) must correspond to the ther- 
mal component only, thus an estimate of the synchrotron 
emission within the beam must be removed from the total 
continuum. 

Fig. [5] shows the distribution of the total continuum, 
the free-free emission estimated from the RRLs with the 
Local To of 7000 K and the difference, that corresponds to 
the synchrotron emission, at fe = 0°. The arrows indicate 
three regions of decrease in the derived synchrotron emis- 
sion coincident in shape and position with Hll regions. This 
is due to an overestimation of the free-free caused by an 
overestimation of To. The opposite occurs for example for 
the Hll region W40, where an increase in the synchrotron 
distribution is correlated with the free-free emission. This in- 
dicates that the electron temperature of this local Hll region 
is higher than 7000 K. Fig.[5]also shows that the synchrotron 
longitude distribution is relatively smooth. We use this fact 
to estimate the synchrotron emission underlying Hll regions 
situated away from any strong synchrotron emitting source. 
The synchrotron background is taken as the average of its 
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Figure 5. Comparison between the total continuum, free-free estimated from the RRLs with To = 7000 K and the synchrotron, versus 
longitude at 6 = 0° . The arrows show three regions of decrease in the derived synchrotron emission coincident in shape and position with 
Hll regions. 



longitude distribution across an Hll region and either side 
of it. 

We selected 16 Hll regions in the area under study that 
have well-defined longitude profiles and spectra with only 
one velocity component, for which To is obtained from equa- 
tion ([31 . Using the longitude profiles taken at the central 
latitude of each Hll region, we estimate the object's longi- 
tude extent from its free-free distribution. The synchrotron 
background, Tsync, is the average of the synchrotron distri- 
bution either side of the Hll region. The thermal continuum 
Tb is calculated from the difference of the total Tc and the 
synchrotron temperatures. These are measured at the centre 
of the object, along with the line integral J Ti^AV in equa- 
tion ([3]). The uncertainties on these quantities are the rms of 
the corresponding maps measured away from the Galactic 
plane, typically at \b\ = 3° —4°. Table[T]lists the Hll regions 
with the corresponding RRL velocities and Galactocentric 
radii, along with the de rived Te . The Galactocentri c distance 
is calculated using the iFich. Blitz fc Starlj (|l989h rotation 
curve, with Ro = 8.5 kpc, assuming that the object moves 
around the Galactic centre in a circular orbit with radial 
velocity V. 

The electron temperatures obtained for the 16 Hll re- 
gions are shown in Fig. [5] as a function of the Galactocentric 
distance. The best linear fit to the data gives: 



To = (3423 ± 472) + (517 ± 99)i?G 



(5) 



where the uncertainties on Te are taken into account by the 
least-square fitting procedure. As shown in Fig. [U] this result 
is in g enera l agreement with th e iPaladini. Davies fc DeZottil 
l|2004h andHh aver et al.l (| 19831 ) results with larger samples. 
The difference between the three lines is mainly caused by 
the lack of data points beyond the solar radius in the present 
sample and the two higher Te values from the local Hll re- 
gions W40 and W45. There is an intrinsic scatter at any 
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9 8 7 6 5 4 3 

Galactocentric radius (kpc) 

Figure 6. The Tc — Rq relationship derived from the results 
in Table [T] The black line gives the best linear fit to the data 
points (equation l[5|l) and the b lue a nd green hn e s are the 
iPaladini. Davies fc DeZottil ||2004|) and IShaver et al] lll983l') re- 
sults, respectively. 



given Rg in all studies due to the different properties of 
each Hll region. 



4.2 The total free-free brightness temperature 

We are now in a position to convert the RRL integral into 
free-free brightness temperature at each pixel using the de- 
rived Tc — Rg relationship in equation ((S)). The Galacto- 
centric distance is calculated using the velocity at peak of 
each spectra. The maximum and minimum of the derived 
Tc map for this region are 4900 K and 8590 K, respectively. 
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Table 1. List of Hll regions selected for the study of the electron temperature v ariation with Galactocentric distance. Rq, in column 
3 is calculated using the RRL central velocity and the lFich, Blitz &: Starld (jToS^) rotation curve with Ro = 8.5 kpc. The velocities and 
corresponding uncertainties, in column 2, are the result of a Gaussian fit. Tq, Tgync, Tff, J Tj^AV are obtained from the longitude 
profiles using the method described in the text; the uncertainties, taken as the rms measured on the maps away from the Galactic plane, 
are 0.6 K, 0.6 K, 0.2 K an d 0.5 K km s~^, resp ectively. To is calculated using equation (|3J. The last column identifies each Hll region by 
its commonly used name jPaladini et alj|2003l) . 



Hll Region 


V 




Tc 


T 

^ sync 


'^ff 




To 


Name 




( km s^i) 


(kpc) 


(K) 


(K) 


(K) 


(K km s-l) 


(K) 




G40.5+2.5 


25.9 + 0.7 


7.1 


9.5 


5.0 


4.5 


1.3 


8360 + 1590 


W45 


G37.8-0.1 


59 + 1 


5.9 


21.1 


6.4 


14.7 


6.0 


6180 + 340 


W47 


G37.7-0.2 


60 + 1 


5.9 


22.3 


7.3 


15.0 


5.8 


6520 + 350 




G36.3-1.7 


62.2 + 0.5 


5.7 


7.0 


3.8 


3.2 


1.3 


6210 + 1480 


RWC179 


G35.6+0.0 


53.4 + 0.3 


6.0 


19.4 


8.8 


10.6 


4.0 


6650 + 410 




G33.1-0.1 


98.2 + 0.4 


4.7 


15.7 


8.6 


7.1 


3.3 


5560 + 530 


Ke78 


G30. 9-0.1 


97.5 + 0.6 


4.6 


61.0 


8.9 


52.1 


23.7 


5660 + 190 


W43 


G30.7-I-0.0 


95.6 + 0.6 


4.6 


77.3 


10.4 


66.9 


30.9 


5580 + 170 


W43 


G30.0-0.1 


97.3 + 0.4 


4.5 


61.0 


8.9 


52.1 


23.7 


5780 + 280 




G28. 8+3.5 


1.5 + 0.4 


8.3 


24.7 


6.0 


18.7 


5.1 


8770 + 560 


W40 


G27.1+0.0 


92.9 + 0.5 


4.4 


20.9 


8.6 


12.3 


4.8 


6500 + 440 


3C387 


G25. 8+0.2 


107.6 + 0.4 


4.0 


27.7 


10.0 


17.6 


7.4 


6060 + 290 




G24.8+0.1 


108.0 + 0.4 


3.9 


28.9 


8.4 


20.5 


9.1 


5790 + 240 


W42 


G23. 5+0.0 


92.2 + 0.4 


4.1 


27.4 


12.0 


15.4 


7.4 


5400 + 280 




G22.9-0.3 


74.8 + 0.4 


4.5 


33.7 


12.5 


21.2 


9.2 


5870 + 220 


W41 


G20. 7-0.1 


53.6 + 0.3 


5.0 


23.2 


10.3 


12.9 


6.3 


5340 + 290 


Ke68 



The average Te at 6 = 0° over the longitude range studied is 
6080 ±970 K. The fact that we are ignoring a possible second 
weaker velocity component does not affect the resulting Tc 
by more than its uncertainty. 

The final free- free map is shown in Fig.[71 along with to- 
tal and synchrotron emission maps for the same region. The 
synchrotron map will be discussed in Section[5l The free-free 
map shows the thermal emission emanating from particular 
objects, mostly lying close to the Galactic plane, whereas the 
diffuse component of the total continuum is mainly associ- 
ated with the non-thermal emission. While there is evidence 
of incomplete separation, e.g. in the W40 Hll region, most of 
the compact sources are well separated as either Hll regions 
or SNRs, and that, along with the diffuse emission, provides 
a visual indication of the reliability of the method. 



4.3 Catalogue of Hll regions 

A catalogue of Hll regions from the present free-free map 
was created using the sourc e extractor program SExtrac- 
tor (|Bertin fc Arnoutd [l99^ . The main parameters to ad- 
just in SExtractor are the background mesh size, detection 
threshold and minimum area which are set to 16, 2a and 
2 pixels, respectively. The filter function used is a Mexi- 
can hat filter of width equal to 3 pixels. The mesh size is 
equivalent to ~ 1° and the width of the filter function is 
similar to the 14.8 arcmin beam. The estimated background 
rms by SExtractor is 0.26 Jy/beam, so the threshold limit is 
0.53Jy/beam. The contrast parameter, which controls the 
deblending of multiple sources, is set to 0.005. These param- 
eters were tuned to recover objects on the Galactic plane 
with sizes consistent with those from higher resolution sur- 
ve ys, such has the 4. 3 arcmin resolution survey at 2.7 GHz 
bv lReich et all (I199G| ). 



Table [5] lists the 57 Hll regions extracted from the free- 
free map. After the source is detected in the filtered map, 
it is fitted for an elliptical shape, using the second order 
moments, whose parameters a and b correspond to the el- 
lipse major and minor axis, respectively. The position angle 
(P.A.) listed is the angle between the semi-major axis of the 
source and the longitude axis, measured counter-clockwise. 
The peak flux of the source is measured directly from the 
input map and subtracted by the background. The uncer- 
tainty on the peak flux is estimated using the rms uncertain- 
ties on the free-free and background maps. Since SExtractor 
does not provide flux densities from Gaussian fitting, we use 
the peak flux and the observed size (a x b) to calculate it 
assuming a Gaussian profile for the source: 

where Sp is the peak fiux in Jy/beam and fwhm is the beam 
FWHM of 14.8 arcmin. The uncertainty on S is estimated 
using the uncertainties on Sp and on a and b, which are given 
by SExtractor. The maximum uncertainty on the observed 
size is 4 arcmin. A source of 5 arcmin in size is observed by 
the 14.8 arcmin beam with a diameter of 15.6 arcmin. This 
small broadening of the beamwidth is difficult to measure, 
especially for a weak object within the noise. For the ob- 
jects whose apparent area is smaller than the beam area, 
the flux density is the same as the peak, S = Sp. In such 
cases, the P.A. has a large uncertainty. The flag parameter 
in Table [5] indicates if: 1 the object has bright neighbours; 2 
the source was originally blended with another one; 3 that 
cases 1 and 2 apply. It shows that deblending occurs mainly 
in the Hll complexes W47, W43 and W42. 
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Table 2. Hll regions in the region I = 44° to 20°, extracted from the free-free map using SExtractor. Column 1 numbers each object; 
columns 2 and 3 are the coordinates; columns 4 and 5 give the angular size of the source, thus a and b deconvolved with the 14.8arcmin 
beam, which is p when the object is less than a few (~ 5) arcmin; column 6 is the angle measured counter-clockwise between the semi- 
major axis of the source and the longitude axis; columns 7 and 8 give the peak and flux density at 1.4 GHz; column 9 gives the internal 
flags generated by SExtractor (see text); column 10 gives the R RL centra l vel ocity from a one or two component Gaussian fit; column 
11 identifies a source by its commonly used name as in the Paladl ni et al.l ll2003i ) catalogue. 



Number 


£ 


b 


6»a 


db 


P.A. 


Peak 


Flux 


Flag 


V 


Name 




(°) 


n 


(') 


(') 


(°) 


(Jy/beam) 


(Jy) 


km 






1 


43.18 


0.00 


10.4 


6.3 


-45 


25.6 ±3.6 


41.2 ± 6.7 





9.4 


W49A 


2 


42.54 


-0.21 


9.7 


P 


-76 


6.0 ± 1.3 


6.8 ± 1.8 





20.4,67.3 




3 


42.06 


-0.58 


1.7 


P 


-9 


3.7 ±0.9 


3.7 ± 0.9 





66.6 




4 


41.63 


0.10 


4.9 


P 


-19 


4.7 ±0.9 


4.7 ± 0.9 





14.9,53.6 




5 


41.08 


-0.19 


21.7 


P 


11 


8.0 ± 1.2 


10.1 ± 2.0 





58.8 




6 


40.54 


2.51 


6.1 


P 


68 


4.5 ± 1.0 


4.5 ± 1.0 





26.2 


W45 


7 


40.00 


-1.37 


P 


P 





4.1 ± 1.2 


4.1 ± 1.2 





47.7 


S74 


8 


39.38 


-0.17 


6.8 


P 


33 


5.0 ± 1.3 


5.0 ± 1.3 


2 


19.4, 60.9 


NRA0591 


9 


38.23 


-0.14 


12.2 


P 


19 


12.4 ± 2.9 


12.4 ±2.9 


1 


56.1,86.9 




10 


37.84 


-0.25 


17.1 


1.7 


47 


17.8 ±3.9 


27.3 ±6.5 


3 


60.7 


W47 


11 


37.60 


-0.01 


32.2 


P 


-43 


16.6 ±3.3 


32.0 ± 7.3 


2 


52.2,86.7 


W47 


12 


37.21 


-0.12 


5.5 


P 


24 


12.0 ±2.5 


12.0 ±2.5 


3 


42.5,80.2 




13 


36.91 


-0.16 


P 


P 


46 


6.2 ± 2.0 


6.2 ±2.0 





40.9,80.4 




14 


36.50 


-0.20 


P 


P 





3.3 ± 1.8 


3.3 ± 1.8 





46.6,73.9 




15 


36.34 


-1.67 


2.4 


P 


-9 


4.3 ± 1.1 


4.3 ± 1.1 





62.1 


RCW179 


16 


35.65 


-0.81 


5.2 


P 


45 


5.0 ± 1.7 


5.0 ± 1.7 


2 


58.8 




17 


35.59 


-0.01 


10.1 


7.7 


86 


10.5 ± 2.4 


14.4 ±3.6 





53.4 




18 


35.27 


1.03 


P 


P 





3.8 ± 1.2 


3.8 ± 1.2 





76.0 




19 


35.19 


-1.69 


27.4 


P 


53 


7.9 ± 1.4 


13.7 ±3.0 





45.3 


W48 


20 


34.88 


0.00 


15.9 


11.6 


69 


9.2 ± 2.5 


17.1 ± 5.0 





53.3 




21 


34.74 


-0.77 


P 


P 


32 


3.2 ± 1.7 


3.2 ± 1.7 





47.3 


W44 


22 


34.52 


-1.05 


P 


P 


47 


3.8 ± 1.2 


3.8 ± 1.2 





48.1 




23 


34.27 


0.13 


19.9 


P 


2 


12.3 ±3.0 


16.5 ± 4.4 


1 


57.2 


W44 


24 


33.57 


-0.06 


21.2 


P 


51 


7.7 ± 2.4 


8.2 ±2.9 





58.9,99.7 




25 


33.14 


-0.10 


16.5 


7.9 


-45 


9.6 ±2.4 


16.4 ± 4.4 





98.3 


Ke78 


26 


32.22 


-0.17 


12.9 


P 


-15 


7.7 ±2.3 


7.7 ± 2.3 


1 


39.9,88.0 




27 


32.17 


0.10 


2.4 


P 


70 


6.6 ±2.0 


6.6 ±2.0 


1 


44.1,90.6 




28 


31.93 


-0.37 


P 


P 


2 


8.4 ±2.5 


8.4 ± 2.5 


1 


42.9,88.7 




29 


31.90 


1.47 


P 


P 





2.2 ± 0.6 


2.2 ±0.6 





56.7 


RCW177 


30 


31.49 


-0.21 


6.1 


P 


-9 


15.1 ±3.6 


15.1 ± 3.6 


3 


39.0,94.7 




31 


30.86 


-0.74 


4.9 


P 


12 


8.4 ±2.8 


8.4 ± 2.8 





94.8 




32 


30.81 


-0.02 


29.1 


12.2 


-11 


82.4 ± 16.3 


236.3 ±50.6 


1 


97.6 


W43,M51,Ke76 


33 


30.27 


-0.19 


31.9 


P 


26 


32.4 ±6.4 


66.0 ± 14.7 


3 


99.8 




34 


29.93 


-0.05 


18.2 


10.6 


83 


24.2 ±4.8 


47.5 ± 10.3 


2 


97.2 




35 


29.21 


0.08 


14.0 


P 


-29 


10.0 ±3.2 


10.1 ± 3.6 





56.9,96.1 




36 


29.09 


-0.65 


P 


P 


-73 


4.1 ± 2.4 


4.1 ± 2.4 





50.0 


RCW175 


37 


28.80 


3.49 


10.9 


5.8 


-52 


19.0 ± 3.0 


30.6 ±5.4 





1.5 


W40 


38 


28.74 


0.17 


24.9 


11.3 


-43 


13.1 ± 3.7 


32.0 ±9.6 


1 


97.2 




39 


28.18 


-0.02 


29.0 


P 


-34 


10.1 ±3.1 


17.9 ±6.0 


1 


94.9 




40 


27.41 


-0.13 


6.8 


P 


-18 


7.8 ±3.0 


7.8 ±3.0 


3 


92.6 


3C387 


41 


27.17 


0.03 


10.6 


P 


-75 


13.0 ±4.5 


10.7 ±3.3 


2 


93.7 


3C387 


42 


26.66 


-0.14 


21.7 


P 


2 


14.9 ±3.9 


22.7 ±4.5 


1 


73.0, 101.7 




43 


26.13 


-0.07 


P 


P 





10.8 ±4.3 


10.8 ±4.3 


3 


103.2 




44 


25.77 


0.14 


25.0 


9.9 


-78 


16.5 ± 5.5 


39.1 ± 13.4 


3 


46.8, 106.7 




45 


25.39 


-0.20 


14.5 


10.6 


-64 


25.0 ± 7.2 


43.1 ± 13.0 


3 


61.4, 109.5 


W42,Ke72 


46 


25.22 


0.31 




n 


18 


8.0 ± 5.8 


8.0 ± 5.8 





44.3, 105.5 


ROW 173 


47 


24.58 


0.17 


33.9 


10.1 


19 


24.0 ± 7.8 


72.5 ± 24.3 


3 


108.0 


W42 


48 


24.55 


-0.16 


29.2 


7.1 


-2 


25.4 ±7.9 


62.1 ± 20.2 


3 


102.2 


W42 


49 


23.62 


0.22 


17.1 


8.6 


53 


19.1 ± 5.7 


33.6 ± 10.5 


3 


90.9 




50 


23.43 


-0.17 


30.6 


5.2 


-30 


22.4 ±5.8 


54.8 ±15.0 


3 


60.2,93.0 


W41 


51 


23.03 


0.60 


P 


P 





2.0 ±2.2 


2.0 ±2.2 





27.7,96.6 


S57 


52 


22.93 


-0.35 


24.7 


15.9 


-13 


23.3 ± 5.3 


66.7 ± 16.3 


3 


75.4 


W41 


53 


22.40 


-0.00 


12.4 


P 


90 


7.6 ±2.2 


7.6 ±2.2 


2 


85.0 




54 


22.22 


-0.17 


7.7 


P 


6 


8.5 ± 2.3 


8.5 ±2.3 


3 


82.1 




55 


21.78 


-0.17 


17.9 


14.2 


37 


8.0 ± 1.8 


17.4 ±4.3 





79.5 


RCW168 


56 


20.74 


-0.10 


13.6 


8.1 


-14 


17.1 ±3.3 


26.5 ± 5.6 





52.8 


Ke68 


57 


20.27 


-0.77 


P 


P 


90 


2.9 ±0.7 


2.9 ±0.7 





48.1 
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Figure 7. Maps of the total continuum, free-free and synciirotron emission at 1.4 GHz and 14.8arcmin resolution. The free-free is 
estimated from the RRL integral using the Tc — Rg relationship from equation JSj. The synchrotron is the difference between the total 
continuum and the free-free emission and shows a narrow diffuse emission confined to the plane. The colour scale is linear and in units 
of brightness temperature (K) . 



4-3.1 Properties of the catalogue 



From Section 13.31 it follows that the RRL integrated fiuxes 
for a given source are accurate to the 10 per cent level. 
However, the ~ 500 K uncertainty on the relationship 
in equation ((5]) means that the integrated fiuxes in Table 
[5] are accurate to the 15 per cent level. As an example, we 
comp are the present result for W40 with the higher reso- 
lution IPaladini et al.l |2003) catalogue, which lists the flux 
density a t 2.7 GHz for 264 Hll regions in this longitude 
range. In iPaladini et al.. the flux density and size of W40 
are 34.0 ± 3.4 Jy and 5.5 ± 2.0 arcmin, respectively. We as- 
sume this Hll region is optically thin at 1.4 GHz based on 



the spectral index a = 0.0, where S oc , given by the 
lAltenhoff etHI (|l97(]| ') flux measurements at 1.4, 2.7 and 
5 GHz with a ~ 10 arcmin beam in each case. Therefore, 
the fiux density at 1.4 GHz extrapolated using the optically 
thin spectral index a = —0.1 is 37.3 ± 3.6 Jy, compared to 
30.6 ± 5.4 from this work. The two results are consistent. 

In a complex of Hll regions the 14.8 arcmin beam cannot 
resolve the individual components, thus comparing the total 
integrated flux with the individual flux densities from the 
IPalad ini ct"aL (200 i s not straightforward. The 2.7 GHz 
survey by Reich et al] ((l990) , at a resolution of 4.3 arcmin, 
was used to identify the individual Hll regions within the 
source G26. 66-0.14 listed in Table [1 These are G26.6-0.1 
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and G26.6-0.3, listed in the IPaladini et al.l catalogue with 
angular sizes of 5.5 ± 1.9 arcmin and 11.0 ± 3.6 arcmin, and 
velocities 104.2 ± 1.8 km s"^ and 69.2 ± 1.6 km s~\ respec- 
tively. The present RRL spectrum at {t, b) = (26? 66, -0?14) 
has two velocity components, 101.7 and 73.0 km s~^, which 
means that the emission is indeed arising from these two 
Hll regions, r. 



4.4 The latitude distribution of individual 
Hll regions 

In this Section we compare the latitu de distribution of the 
Hll regions from the higher resolution IPaladini et al.l (|2003l ) 
catalogue with that from the prese nt, lower reso lution cat- 
alogue. There are 159 objects in P aladini et al.l with RRL 
velocity measurements, in this longitude range. The listed 
2.7 GHz flux densities are extrapolated to 1.4 GHz assuming 
the free-free is optically thin therefore well approximated by 
a power law with a = —0.1. The flux densities and sizes 
from both catalogues are used to simulate Gaussian profiles 
and to create the spatial distribution of the Hll regions in 
the form of brightness temperature maps at 1.4 GHz and 
14.8 arcmin resolution. We have also estimated the contri- 
bution of UCHll regions to the total free-free emission in 
the area under study using the apparent brightness temper- 
ature of the 220 catalogued compact Hll r egion s at 1.4 GHz 
and 5 arcsec resolution by iGiveon et all l|2005l ). About 20 
per cent of the objects are likely to be optically thick at this 
frequency. However the contribution from the UCHll regions 
is found to be relatively low, 1 — 2 per cent. Thus the op- 
tically thin assumption for the free-free emission is not sig- 
nificantly affected by the presence of these compact objects. 
The UC objects contribute with 1 — 2 per cent to the total 
free-free emission, therefore they is not likely to affect the 
With Therefore the contribution of UC objects is small to 
affect the optically thin assumption on the spectral index 
of free-free emission. The comparison is shown in Fig. |8] (a) 
as a function of latitude, averaged over the whole longitude 
range. The two distribu tions are very similar, with FWHM 
of 0?48 and 0?53 for the IPaladini" et al.l and the present lists, 
respectively. The catalogue from this work includes the more 
extended emission around and between fe atures which is not 
included in the higher resolution Paladini et al.l l|2003l ) cat- 
alogue, which is missing ~ 20% of the flux. 

Fig. (Sj (b) compares the latitude distribution of the 
Hll regions recovered with SExtractor with that of the to- 
tal free-free emission, in the longitude range 20° to 44°. It 
can be seen that the Hll regions are a narrower distribution 
than the total diffuse emission, with about half the FWHM. 
The individual Hll regions account for ~ 30% of the total 
free-free emission in this region of the Galaxy. This result 
is similar to that found in Paper I for the smaller longitude 
range £ = 36° to 44°. 



20°< ? <44° 
18 1 ' ■ ' ' ' ' ^VT' ' 




3 I ^ . ^ I ^ ^ ^ I ^ ^ ^ \ ^ ^ . I 

-4 -2 2 4 

Galactic Latitude (°) 

Figure 9. The synchrotron emission versus latitude, averaged 
over the 24° ^-range using 4 arcmin latitude bands, obtained with 
the Tc - Rg relationship (full), Tc = 7000 K (dotted) and Tc = 
5000 K (dashed). 

The result in show in Fig. [7] Well-known SNRs are iden- 
tified in the map along with a broad background emission 
falling to ~ 4K at |&| ;$ 4°. There are some positive and 
negative residuals at the position of a few Hll regions in the 
synchrotron map. W40 and W48 are likely to be underes- 
timated in free-free emission caused by the adoption of an 
electron temperature for the Hll region less than the value 
used for the diffuse emission. In the case of W40, the broad 
emission underlying the compact source is thought to be of 
non-thermal origin. The opposite applies to W47 and W43. 
These are groups of several Hll regions, unlikely to be char- 
acterised by a single To value. 

Fig. [9] gives the synchrotron emission versus latitude, 
averaged over the 24° range in longitude, using 4 arcmin 
latitude bands. The three lines are the different synchrotron 
results using To = T^Rq), Te = 7000 K and To = 5000 K 
in the free-free estimation. The similarity between the three 
plots for l&l 1° reflects the fact that the RRL line integral 
is low away from the plane, and therefore the synchrotron 
emission is the major contributor at 1.4 GHz for |6| 1°. 
The contribution of the strongest SNR in this region, W44 
(G34. 7-0.4), is visible at b ^ -0°4. The smoother distri- 
bution given by the full and dashed lines in Fig. [5] with 
relation to the dotted line, seems to indicate that Te of the 
diffuse ionised gas on the Galactic plane is similar to that of 
the individual Hll regions and not higher. The synchrotron 
distribution appears to have two components - one with a 
FWHM of 10° and a narrow Galactic plane component 
with a FWHM of 2° . This narrow component is identified 
here for the first time after direct subtraction of the free-free 
emission and is presumably arising from the SNRs on the 
plane which have different ages. 



5 DERIVATION OF A 1.4 GHZ 
SYNCHROTRON TEMPLATE 

The synchrotron emission at 1.4 GHz for the region £ — 
20° to 44° and |&| < 4° can now be obtained by subtract- 
ing the final free-free emission from the total continuum. 



5.1 Catalogue of synchrotron features 

The synchrotron map derived in this work, as the best mea- 
sure of synchrotron emission at 1.4 GHz, enables the de- 
tection of synchrotron features, most of which are SNRs. 
SExtractor was used to create a list of SNRs in the region 
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Galactic latitude {°) Galactic latitude {°) 

Figure 8. (a) Comparison between the contribution from the 159 iPaladini et ahl l l2003h sources (dotted line) and the 57 sources from 
extracted from this work (full line), versus latitude, (b) Separation between the diffuse emission and the emission from individual 
Hll regions at 1.4 GHz. The full line represents the total free-free emission and the dotted is the contribution from the SExtractor 
Hll regions. Both plots are for the longitude range 20° to 44° with 0?2 latitude bins. 



I = 44° to 20°. The background mesh size, the filter func- 
tion and the detection minimum area are the same as for 
the free-free map but the detection threshold is decreased 
from 2 to l.Scr. The rms of the resulting background map 
is 0.69 Jy/beam and the threshold limit for a detection is 
therefore 1.03 Jy/beam. A Gaussian filter was also used in 
order to detect diffuse and extended such as the W50 SNR, 
not easily detected by the Mexican hat filter. Table [3] gives 
the angular size, position angle, peak and total fiux density 
for the 26 SNRs recovered by SExtractor. 

The latest publish ed compilation of Galactic SNRs is 
given bv lGreenI (|2009al '). It contains 36 SNRs in the ^-range 
under study, 23 of which are detected in this surv ey and 
listed in Table [31 Objects 3, 4 and 10, are not in the iGreenI 
catalogue and are described in mo re d etail in the next sec- 
tion. Most of the 13 SNRs in the [Green catalogue not de- 
tected in this work have flux densities below 2.7 Jy, which is 
comparable with the lowest flux density in Table[3l Moreover 
some of the SNRs are poorly defined, without clear angular 
sizes or fiux measurements. 

Most of the sources in Table [3] are smaller than the 
beam, their sizes and fi ux dens ities are in overall agreement 
with the values in the iGreenI catalogue. For example the 
strongest SNR in this region of the Galaxy, W44, is fitted in 
the present synchrotron map with a size of 26'.0 x 15'. 4 and a 
fiux de nsity of 208 ± 19 Jy. The fiux density given bv lGreenI 
(|2009al 'l for this object is 203 Jy at 1.4 GHz, extrapolated 
from 1 GHz with the listed spectral index a = 0.37, and the 
angular size is 35' x 27'. 



5.2 Comments on individual synchrotron sources 

5.2.1 G42.12-0.21 and G41. 45+0.39 

These ob jects correspond to two of the three SNR candi- 
dates by Ikaplan et all l|2002l ). The angular sizes and fiux 
densities found here are larger than the values obtained by 
iKaplan et all (|2002|) using the 1.4 GHz NRAO VLA Sky Sur- 
vev fNVSS. ICondon et al] Il998l ) data, which may resolve 



out some of the emission. This can also be due to the fact 
that these two faint objects are recovered with the Gaus- 
sian filter, thus are likely to be confused with the strong 
background on the plane. Moreover, they are both fiagged 
as having bright sou rces nearby. The third possible SNR in 
iKaplan et aTT(|2002| ), G43.5-h0.6, is too faint in the NVSS to 
derive a flux density and it is indeed only marginally visible, 
at a ~ 1.2a level, in the present synchrotron map. 



5.2.2 G 3 5. 6- 0.4 



This SNR has been recently re- identified b vlGreenI ll2009b ) 
using VLA Galactic Plane Survey (VGPS. Istil et al.ll2006l ) 
data and data from other single-dish surveys at higher 
frequency to confirm the non-thermal spectral index of this 
object. G35.6-0.4 is seen in the VGPS map with a size of 
~ 15' X 11' and has an fiux density of 7.8 Jy. These results 
are consistent with what is found in this work, 6.7 ± 1.5 Jy 
for a relatively compact source of size less than or equal to 
that of the beam. 



We have thus confirmed the synchrotron nature of the 
objects G42. 12-0.21, G41.45-f0.39 and G35. 6-0.4 using the 
present survey. 

The simila r resol ution of the 2.7 GHz survey 
by [R eich et al' ('l99Cj)_ and the 100 /^m data by 
iMivill e-Deschcncs & La gachd (|2005) allows the in- 
vestigation of the spatial coincidence between dust 
and radio continuum emission, used to distinguish be- 
tween Hii regions and SNRs jHaslam fc Osbornd Il987l : 
iBroadbent. Osborne fc HaslamI 1 19891 ). This correspondence 
helps us understand that the artifacts seen in the syn- 
chrotron map around the Hii region complexes W47, W43 
and W42 are correlated with individual Hll regions, and 
are thus caused by electron temperature variations. On the 
other hand, the lack of significant dust emission for the 
sources at (^,6) = (30?1,-H?3) and {l,h) = (21?0,-f2?0) 
in the synchrotron map, points towards a non-thermal 
origin. Both sources are extragalactic, even though only 
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Table 3. SNRs in the region £ = 44° to 20°, extracted from the present synchrotron map using SExtr actor. The colu mn descriptions 
are the same as in Table [2] The last column identifies a source by its commonly used name as in the iGreenl ll2009al l catalogue. The 
Flag parameter 16 for the first and last SNR listed means that the objects are close to the image boundary. The 7 objects with the star 
symbol before the number are those extracted using a Gaussian filter instead of the Mexican hat filter. 



Number 


i 


b 


da 


6b 


P.A. 


Peak 


Flux 


Flag 


Name 




(°) 


n 


(') 


(') 


n 


(Jy/beam) 


(Jy) 






*1 


43.98 


1.60 


35.4 


25.4 


-77 


2.0 ± 1.1 


10.1 ± 5.8 


16.0 




2 


43.23 


-0.13 


16.0 


V 


58 


29.8 ± 1.3 


49.7 ± 5.5 


0.0 


W49B 


*3 


42.12 


-0.21 


51.34 


41.1 


78 


3.8 ± 1.3 


40.3 ± 14.4 


1.0 




*4 


41.45 


0.39 


26.8 


20.3 


1 


4.3 ± 1.4 


15.0 ± 5.1 


1.0 




5 


41.10 


-0.30 


P 


P 


-3 


13.3 ± 1.5 


16.0 ± 1.8 


0.0 


3C397 


*6 


40.58 


-0.45 


34.0 


28.0 


-13 


5.0 ± 1.6 


26.9 ±9.2 


3.0 




*7 


39.78 


-2.34 


113.4 


53.2 


73 


4.9 ± 1.3 


141.3 ± 40.3 


17.0 


W50 


8 


39.22 


-0.32 


4.2 


P 


-14 


12.7 ± 1.8 


12.7 ± 1.8 


0.0 


3C396 


*9 


36.61 


-0.83 


21.0 


P 


-51 


2.3 ±1.4 


4.0 ± 2.4 


0.0 




10 


35.59 


-0.44 


P 


P 


55 


6.7± 1.5 


6.7 ± 1.5 


0.0 




11 


34.68 


-0.44 


26.0 


15.4 


-56 


71.2 ± 1.6 


208.1 ± 18.6 


0.0 


W44 


12 


33.67 


0.06 


P 


P 


43 


10.3 ± 1.6 


10.3 ± 1.6 


0.0 


Kes79 


*13 


33.21 


-0.51 


7.9 


P 


14 


2.3 ± 1.4 


2.3 ± 1.4 


0.0 




14 


32.86 


-0.05 


P 


P 


-1 


6.5 ± 1.5 


6.5 ± 1.5 


0.0 


Kes78 


15 


31.89 


0.04 


3.9 


P 


25 


14.7 ± 1.7 


14.7 ± 1.7 


0.0 


3C391 


16 


29.69 


-0.27 


P 


P 


80 


7.9 ± 2.0 


9.5 ± 2.4 


0.0 


Kcs75 


17 


28.61 


-0.05 


P 


P 


-72 


8.6 ±1.8 


8.6 ± 1.8 


0.0 




18 


27.75 


0.60 


P 


P 





10.9 ± 1.4 


10.9 ± 1.4 


0.0 




19 


27.37 


0.00 


P 


P 


-51 


8.9 ± 1.3 


10.7 ± 1.6 


3.0 


4C-04.71 


20 


24.73 


-0.70 


P 


P 


90 


5.4 ±1.5 


5.4 ± 1.5 


0.0 




21 


24.58 


0.62 


10.1 


P 


23 


11.4 ± 1.7 


12.7 ±3.2 


0.0 




22 


23.17 


-0.24 


8.8 


P 


81 


13.0 ±2.3 


14.5 ±4.1 


2.0 


W41 


23 


22.67 


-0.24 


17.3 


2.4 


73 


15.5 ± 2.0 


24.1 ± 5.1 


0.0 




24 


21.81 


-0.57 


14.5 


7.3 


34 


32.5 ± 2.0 


51.0 ±6.2 


0.0 


Kes69 


25 


21.50 


-0.87 


P 


P 





8.2±2.1 


9.8 ± 2.5 


0.0 




26 


19.94 


-0.21 


17.1 


P 


-47 


12.0 ± 1.7 


16.4 ±4.3 


16.0 






Figure 10. (a) Contour map of the radio source {I, b) = (30? 13, +1?34) from the VGPS survey at 1.4 GHz and 1 arcmin resolution. The 
contours are given at every 5% from 5 to 90% of 0.92 Jy/beam, which corresponds to the maximum on the southern lobe. The first 4 
contours are labelled, (b) The Hi spectrum from the VGPS survey against the radio source, taken at {£, b) = (30?15, 1?31) (full line). The 
dotted and dashed lines are the off source spectra, measured at the same b and £ = 31? 13 and 31? 17, respectively. The thick long-dashed 
line is the absorption spectrum of the source which shows that absorption is detected at all velocities, inside and outside the Solar circle, 
confirming its extragalactic origin. 
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(I, b) = (2 1?0, +2?0) is currently identified as such in the 
hterature (|Clark fc Crawford! 1 197^ ). 



5.2.3 The double-lobed radio source (£, b) = (30?13, +1?34) 

The source at (£, b) = (30? 13, +1?34) has been detected in 
previous survey s (IReich et al ][l99 0; Reic h. Reich fc FuerstI 
Il990l : [sie3ll995l : ICohen et al.fl2007i ) but not identified. We 
use VGPS Hi data to identify this object a radio galaxy. The 
contour map of Fig. [TU] (a) is the VGPS 1.4 GHz emission 
of the radio source. At this survey of resolution 1 arcmin 
a double-peaked feature is clearly visible. The elongated 
structure has been detected at lower frequencies, in par- 
ticular at 74 MHz in the VL A Low-Frequency Sky Sur- 
vey (VLSS. ICohen et alllioOTh . The VLSS catalogue gives 
a flux density of 41.97 ± 4.24 Jy and an angular size of 
2'.05 X 0'.59 with a 1.3arcmi r i bea m. At 1.4 GHz, the sur- 
vey bv lReich. Reich fc Fuers^ (|l990l ) at 9.4 arcmin resolution 
gives a flux density of 4.62 ± 0.05 Jy for a slightly extended 
source of angular size les s than 11 arcmin. The higher reso- 
lution survey at 2.7 GHz (jReich et al.lll990l ) detects a source 
of 7'.4 X 4'.3 size and flux density of 3.57 ± 0.04 Jy. Combin- 
ing the three measurements results in a spectral index of 
a = —0 .71 ± 0.05, w hich is within the observed values for 
a SNR l|Green|[2009al ). Nevertheless, the shape of the radio 
source and its angular size points toward an extragalactic 
origm. iHelfand et all (|l992i ). on their 20 cm study with the 
VLA and identification of 1457 compact sources near the 
Galactic plane, also suggest that this is an extragalactic ob- 
ject. If this is the case, such a bright and extended object, 
has probably been missed from extragalactic source cata- 
logues due to its proximity to the Galactic plane. 

The latitude coverage of the VGPS Hi spectral cube, 
|6| J; 1?32, only includes the southern lobe of the source, 
which is enough to obtain and compare spectra on and off 
the source. These are shown in Fig. [10] (b). Hi absorption 
against the object is detected at all velocities including nega- 
tive velocities, which gives the confirmation that it is outside 
the Galaxy. 

This interesting object deserves further study at higher 
angular resolution. 



6 COMPARISON WITH THE WMAP 
MAXIMUM ENTROPY MODEL 

In this Section we extend the analysis performed in Pa- 
per I using the three times larger longitude coverage of 
the present RRL map to compare with the WMAP pre- 
diction of the free-free emission at 23 GHfl The Maxi- 
mum Entropy Model (MEM. [Gold et al.|[201ll ') is a spatial 
and spectral fit that uses external templates in regions of 
low signal-to-noise. The spectral indices for the free-free 
and thermal dust emission are fixed, with /3 = —2.14 and 
/3 = -1-2.0, respectively (Tb oc u^). Any additional compo- 
nent of emission such as the anomalous microwave emission 
is included in the synchrotron component. To have a re- 
liable prediction of the free-free emission on the Galactic 
plane is important not only for foregrounds and component 
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Figure 11. Free-free emission estimated from the RRLs (dashed 
line) and from the WMAP 23 GHz MEM (full line) versus lati- 
tude, averaged over the 24° longitude range. 



separation but also for Galactic st udies such as 3-D Galaxy 
llPlanck Collaboration et al.ll20lil ) and magnetic field stud- 
ies (|jaffe et al.ll201lf ). 

In Paper I the comparison between the free-free esti- 
mated from the RRLs and from the WMAP seven year 
MEM suggested that the electron temperature of the ionised 
medium in the ^-range 36° to 39° is about 8000 K. This 
is based on the higher MEM distribution on the Galac- 
tic plane, by 30 per cent, with relation to the predic- 
tion from the RRLs. Fig. [TT] shows the comparison be- 
tween the two datasets using the extended longitude range 
£ = 20° to 44° from the present survey, at 23 GHz. The 
spectral index used to extrapolate the free-free bright- 
ness temperature from the RRLs is /3 = —2.13 which is 
the average value across the range of electron tempera- 
tures given by equation ([5]). The free- free spectral index 
P is a slow function of frequency an d electron tempera- 
ture (|Dickinson. Davies fc David l2003l ). decreasing as the 
frequency increases. The steeper spectral index used in the 
MEM analysis reflects this behaviour. The RRL data are 
smoothed to the MEM resolution of 1°. The FWHM of 
the MEM latitude distribution is 1?48, slightly broader then 
that of the RRLs, 1?36, with a ~ 50 per cent higher peak. 
The free-free brightness temperature depends on T^'^^ (Sec- 
tion[l|. Therefore this difference implies an electron temper- 
ature of ~ 10000 K for the ionised gas on the Galactic plane, 
compared with the mean 6080 K used in the RRL conversion 
to free-free (Section I4.2|l . Based on the results presented in 
previous sections, we expect the electron temperature of the 
diffuse to be around 6000 K in this region of the Galactic 
plane. 

Similar results are reported bv lJaffe et al.l(l201ir) . They 
recover 80 per cent or less of the MEM free-free emis- 
sion in their model of the Galactic magnetic field. Re- 
sults from the Planck satellite also indicate that the MEM 
overpredicts the ionised compone nt on the Galactic plane 
(jPlanck Collaboration et al.l l201lh. In this case, the MEM 
estimation is ~ 20 per cent higher than that recovered from 
the Galactic inversion model performed. 

The fact that the MEM model does not fit for a sep- 
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Figure 12. The longitude-velocity distribution of RRL bright- 
ness temperature within 1° of the Galactic plane. Contours are 
at —8, —6, —4, —2, —1, 1, 2, 4, 6, 8 and 10 per cent, then at every 
5 per cent from 10 to 50 per cent and every 10 per cent from 60 to 
90 per cent of the maximum line temperature observed, 272 mK. 
The first 6 contours are labelled and the dotted lines correspond 
to the negative contours. 



arate anomalous microwave emission (AME) template may 
contribute to the excess deduced for the free-free template 
relative to the RRL observations. Furthermore, we note that 
the use of a flatter spectral index for the dust emission, with 
respect to fi = -1-2.0, would allow a lower free-free contribu- 
tion at lower frequencies. 

The ratio between the MEM and RRL predictions of 
the free-free emission varies with longitude. The relative con- 
tribution of these two diagnostics will be investigated in a 
future paper, using the full ZOA Galactic plane coverage. 



7 THE 3-D GALAXY 

The present RRL survey is unique for it provides information 
on the distribution of individual Hll regions and the diffuse 
ionised gas on the Galactic plane. 



T.l The 



V diagram 



Fig. ll2l shows the contours of RRL temperature in longitude- 
velocity space, for spectra averaged in 0?5 in longitude 
within 1° of the Galactic plane. The first positive contour 
(full line) is at ~ Icr. It shows that there is no H RRL emis- 
sion at negative velocities in this longitude range: the two 
vertical bands of emission a± V < Q km s~^ are associated 
with the He/C RRLs fSection l3.3p . This res ult is very s imilar 
to that found in the H166a line survey bv lLockmanI (|l976l ) 
at a resolution of 21 arcmin, which detected emission in all 
the 1° intervals from 4° to 44°. Both l — V maps show a min- 
imum of emission at £ ~ 40° , no emission at km s~^ and a 
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Figure 13. Radial distribution of the ionised emission (full line) 
in the region i = 20° to 44° and |b| < 1°. The Galactic dis- 
tance is calculated for the 1503 discrete RRL components using 
-Ro = 8.5 kpc. The dotted line is the ra dial distribution of th e 471 
Hn regions from the co mbination of the lPaladini et al] ||2003| ) and 
lAnderson et al.l l|201lh lists. 



large change in emission with longitude. The main band of 
emission from the H RRLs shows evidence o f a slope in lon- 
gitude, also seen in the Yll l — V diagram in lLockman|[l976l . 
This is the result of circular rotation around the Galactic 
centre, by which gas at higher longitude is further away from 
the centre and thus has lower radial velocity. 



7.2 The radial distribution of the ionised gas 

The RRL velocity, when converted into a distance to the 
Galactic centre can be used to obtain the radial distribu- 
tion of the ionised gas. This is shown in Fig. [TS] The RRL 
spectra are averaged into a 12 arcmin pixel grid, to reduce 
the correlation between the pixels, for the longitude range 
I = 20° to 44° and |6| < 1°. The Gal actocentric radius is 
calculated as in Section 14.11 using the iFich. Blitz fc StarkI 
( 1989) rotation curve, with Rq — 8.5 kpc. The central veloc- 
ity of the RRLs is found by fitting a one or two component 
Gaussian function to the spectra, in the case of a single 
or double velocity component, respectively. From 1331 lines 
of sight we detect 1503 discrete RRL components. Fig. [13] 
also shows the distribution of the 471 Ifll regions in this 
longit ude r ange, from the c o mbin ation of the lPaladini et ahl 
(j2003l) and lAnderson et al.l l|201ll ) samples (scaled by a fac- 
tor of 3). The Gre en Bank Telescop e Hll Region Discov- 
ery Survey (HRDS. lBania et al.ll2010f ) has detected 313 new 
Hll regions between i = 20° and 44°, adding to the 145 
objects with R RL measurements in present region from the 
IPaladini et al.l list. 

The two distributions are very similar and show nar- 
row peaks, « 0.75 kpc wide, at ~ 4.4 kpc and 5.9 kpc that 
correspond to the Scutum and Sagittarius spiral arms, re- 
spectively. This confirms the spatial correlation between 
the diffuse gas and the individual Hll regions, observed by 
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(a) 80 < V < 130 kms/s 
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Figure 14. Maps of RRL emission integrated over 50 km s~^, centred at 105 km (a), 55 km s^^ (b) and 5 km s~^ (c), at a resolution 
of 14.8arcmin. Contours are given at every O.IK km s~^ from —0.4 to 0.4K km s~^, every 0.2K km s^^ from 0.4 to 2 K km s^^ and 
then at 2.5, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25 K km s"!. The negative contours are dotted. 



iHart fc Pediail \l<d7^ ) and lLockmanl \l%7^ ) in the H166a line 
using pointed observations along the Galactic plane. The Rg 
values of the peaks correspond to tangent points in longitude 
of 30° and 45° and terminal velocities of 110 km s~^ and 
63 km s~^, respectively. The fact that the present survey 
covers I < 44° suggests that it is missing part of the emission 
from the tangent point of the Sagittarius arm at ^ = 50° . 

7.3 The z-distribution of the ionised gas 

Having identified the Sagittarius and Scutum arms in the 
Rg and velocity space, we are now able to separate the 
emission from each of the spiral arms and investigate their 
azimuthal extent. Figs. [14] (a) and (b) show the integrated 
RRL emission over 50 km s~^, from 80 to 130 km s~^ and 
from 30 to 80 km s~^, respectively. Fig. [TJ] (a) shows 
the narrow distribution of the high-velocity emission from 
the Scutum arm that covers mainly £ < 34°. The lower 
emission for £ > 34° originates in the wings of RRLs at 
V ^ 80 km s~^. The emission from the Sagittarius arm cov- 
ers the whole longitude range and has a broader distribution 



around the plane. Fig. [14] (c) shows that the RRL emission 
at V ^ km s^-^ is attributed to the W49 and W40 Hll re- 
gions, with some emission from W45, at V ~ 25 km s~^ and 
from the positive latitude extent of the W47 complex. 

The latitude profiles of Figs. [15] (a) and (b) show the 
narrower z-distribution of the Scutum arm, with a FWHM 
of ~ 0?78, compared with the slightly broader distribution 
of the Sagittarius arm, with FWHM of ~ 0?85. The Galacto- 
centric distances to the tangent points give a mean distance 
to the Sun of 6.0 and 7.3 kpc for the Sagittarius and Scutum 
arms, respectively. Using the observed FWHM of the lati- 
tude profiles gives a z-thickness of about 100 pc, a ^ 45 pc, 
for the diffuse ionised gas about the Galactic pl anjf|. This 
compa res with a ~ 70 pc, the value found by iRevnoldsl 
(! l99ll ) using pulsar Dispersion Measures and reflects the fact 
that emission measures trace the densest parts of the ionised 
gas, EM oc rig, compared to dispersion measures which trace 
rio. The latitude distribution of the Sagittarius arm, Fig. 1151 



^ Where a = FWHM/2^21n(2) 
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Figure 15. The latitude profiles that correspond to maps from Figs. [14] (a) and (b), respectively, averaged over the longitude ranges 
indicated. 



(b), shows a broader component of FWHM ~ 2° that corre- 
sponds to the emission from the near side, whereas emission 
from the far side is responsible for the narrow component. 



8 CONCLUSIONS 

This work has confirmed that the HIPASS/ZOA survey can 
provide a reliable estimate of the RRL emission on the 
Galactic plane at 1.4 GHz and 14.8arcmin resolution. The 
use of the total continuum emission from the same survey 
has enabled the electron temperature of individual Hll re- 
gions in the area under study to be estimated, showing the 
well-known Galactocentric gradient consistent with previous 
results. 

The unambiguous determination of the free-free emis- 
sion from the present work is used to derive the synchrotron 
map at 1.4 GHz. The latitude distribution of the synchrotron 
emission shows a narrow component on the Galactic plane 
which is found here for the first time clear of the free-free 
emission by direct subtraction. 

The free-free and synchrotron maps are used to cre- 
ate a list of 57 Hll regions and 26 SNRs, respectively. The 
fact that we d etect a recently re-identified SNR, G35.6-0.4 
l|Greenll2009bl ) illustrates the potential of this survey to de- 
tect and identify new objects. The contribution of the in- 
dividual Hll regions to the total free-free emission is found 
to be around 30 per cent. Their latitude distribution has a 
FWHM of 0?5, about half the width of the total free-free 
distribution. This can be converted into a z-thickness for 
the Hll emission at known distances and is consistent with 
previous results on the distribution of the OB stars and the 
width of the thin disk from Dispersion Measures. 

The present estimate, based on the integrated RRLs 
converted to brightness temperature of the free- free emission 
recovers about 50 per cent of the WMAP MEM prediction. 
Possible reasons for this discrepancy are likely associated 
with the thermal dust spectral index and with the fact that 
there is no separate AME template from MEM. 

This study will soon be extended to the inner Galaxy 



region, £ — —50° to 50°, using the improved data analy- 
sis pipeline described in Section [3.11 We will then focus on 
the different applications of this result, namely on the inves- 
tigation of the spectral energy distribution of the different 
emission components on the Galactic plane and especially 
on the AME, as a function of longitude. Also, the combi- 
nation of Hq and RRL data to complement the free-free 
templates on the Galactic plane, which is the ultimate goal 
of this project, will be subject of a more detailed study. 
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